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The pseudoscalar and vector heavy meson masses are calculated in the e-regime of Heavy Meson 
Chiral Perturbation Theory to order e*. The results presented will allow the determination of low- 
energy coefficients (LECs) directly from Lattice QCD calculations of the heavy mesons masses for 
lattices that satisfy the e-regime criteria. In particular, the LECs that parametrize the NLO volume 
dependance of the heavy meson masses are necessary for evaluating the light pseudoscalar meson 
(tt, K, 77) and heavy meson {{D^\ , Dt} , {B~ , B*^ , B^l}) scattering phase shifts. 

INTRODUCTION 

Understanding the properties of systems composed of heavy mesons, containing a single heavy quark, and the 
pseudo-Goldstone bosons (pGB) of Quantum Chromodynamics (QCD) is currently a topic of high interest. This 
interest has been partly triggered by the renaissance of charmonium and open-charm studies. A resonance that has 
initiated much discussion is the narrow D*g(2317), first observed by the BaBar collaboration yy . This resonance couples 
to the S-wave DK continuum scattering state f2^,'3' . At low energies the strength of the DK interaction is predominantly 
parametrized by the scattering length. This has resulted in several theoretical studies that have attempted to determine 
the S-wave scattering lengths in the pGB-heavy meson scattering channels [1H7] . The determination of these scattering 
lengths would not only help discern the heavy meson spectrum, but is needed in order to evaluate transport coefficients 
of systems containing heavy-light mesonic species, e.g. the hadronic phase of heavy ion collisions. 

Currently, a combination of effective field theories (EFTs) and Lattice Quantum Chromodynamics (LQCD) provides 
the best option for performing reliable calculations of low energy QCD observables (reviews on these topics include 
PfM]). Heavy Meson Chiral Perturbation Theory (HMxPT) [TSHiT] is the low-energy EFT for studying strong- 
interaction quantities of mesons containing a single heavy quark and a single light antiquark. The non-perturbative 
QCD contributions to HM^PT are parametrized by low-energy coefficients (LECs). The predictive power of HMxPT 
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is currently limited by the poor determination of these LECs, e.g. currently next-to-leading order (NLO) LECs are 
determined within a factor of three of precision, resulting in scattering lengths that are known within a factor of 
three [BJ [7] . The results outlined in the work will help reduce the uncertainties of LECs needed in the evaluation of 
pGB-heavy meson scattering. 

Historically, LQCD calculations have used moderate volumes and unphysically large pion masses. With advances in 
computing technology, performing LQCD calculation at the physical point {itIt^ « 140 MeV) of QCD is now a reality. 
However, limited computer resources require state of the art calculation with physical pion masses to be performed 
with small physical volumes. This leads to sizable volume effects contributing to the quantities of interest, and while it 
is natural to want to remove them, these effects can hold physically important information. More specifically, volume 
effects are parametrized by the LECs of the EFT, therefore by evaluating physical observables in small volumes one 
can determine the LECs. 

In an infinite volume, the expansion parameters of HM^PT are p/A^, m;/A^, and Aqcd/"t^q, where p is the 
characteristic momentum of the interaction, m; is the mass of the light Goldstone bosons, mq is the heavy quark 
mass, is the chiral symmetry breaking scale, and Aqcd is the characteristic scale of QCD. In a finite volume this 
expansion scheme is consistent in the p-regime |18l I19j . 

However, for volumes smaller than the Compton wavelength of the Goldstone bosons, the zero momentum mode is 
enhanced with respect to the non-zero modes, and an alternative expansion scheme must be utilized [20]. The regime 
where the pion zero-modes must be integrated over explicitly while still treating the non-zero modes perturbatively is 
known as the e-regime |20H28j . In the e-regime, a new expansion parameter is introduced, e ^ ^ 27r//3A^ and 

TO//A^, where L and /3 are the spatial and temporal extents, respectively. At leading order, one may associate 
Aq, the hyperfine splitting between the pseudoscalar meson P D^, D^}, {B^ , B^}) and its respective vector 

meson P* {{D*° , D*+ , 0'*+}, {B*- , B*° , Bf}), with the physical values on the order of 140 MeV and 50 MeV for 
the charm and bottom mesons respectively. Therefore, it is reasonable to expect the hyperfine splitting to contribute 
at order for charmed mesons (e^ ^ Ac/A^) and approximately at order for bottom mesons {e^ ^ Ab/A^). For 
the sake of generality, both scenarios are considered. 

This study presents the volume dependence of the heavy meson masses at next-to-leading order (NLO), ©(e'*), in 
the "mixed regime" of SU(2) and SU(3) HMxPT. In the "mixed regime", the physical pion mass is small compared 
to the IR cutoff and therefore fall within the e-regime, while the kaon and eta still satisfy the p-regime criteria 



[2Sl[3n]- Therefore in the mixed regime, the expansion in the tt and {K,ri} masses is treated separately in order 
to satisfy ~ m^r/A^ and e '-^ to^/A^ ^ m^/A^. The 0{e^) volume dependence of the heavy meson mass for an 
SU(2)i X SU(2)fl chiral theory with static heavy quarks has been previously calculated [3D]. Unfortunately, the 0{e^) 
result does not give the predictive power necessary to extract the LECs desired. 

HEAVY MESON CHIRAL PERTURBATION THEORY 



The field multiplet of the pseudoscalar P and the vector P* can be conveniently represented as a single field operator 



[HHm, 



1+ ^ 



(1) 



where is the velocity of the heavy meson. The representation of the composite field Ha{x) assures it transforms 
as an SU(2) spinor under heavy-quark spin rotations, under the unbroken SU(3)v' symmetry it transforms as an 
element of the 3 fundamental representation (as denoted by the subscript "a"), and under both Lorentz and parity 
transformations it is a bilinear. In the rest frame, the LO HMxPT Lagrangian in mg and A^ consistent with 
spontaneously broken SU(3)l x SU(3)ii is PW7] : 



£" = -iTr [HaiDyaHb] - g Tr 



HaHb^l ■ Abalb 



f2 

Tr [9^E^9^S] + ^—Tr [TWS^ + h.c] , (2) 



where A4 = ^diag{rri^, m^, 2m\ — m'^) is the light meson mass matrix, 7 is the spatial component of 7'', = d^+V^ 
is the covariant derivative, / is the pion decay constant, and the Goldstone bosons are encapsulated in the operators. 
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S = C = cxp 
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At NLO in HM^PT there are a large number of corrections to the Lagrangian that are consistent with velocity 
reparametrization invariance (VRI) [31) , but the terms that will contribute to the volume dependance of the mass are 
the following: 



HaHbJ ■ Abal5 
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X 



Tr [Haa^^'^Haa^,] + ^Tr {A" A 
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Tr [HaHb i^M^ + /i.c.),J + ^Tr [HaHa {^M^ + h.c\^] + ^Tr [HaH^AlbA^ + j^Tr [HaH.Acb ■ Aba] 
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At leading order, the hyperfine splitting can be written in terms of the LEG A, Aq = 



ZERO-MODES INTEGRATION IN THE e-REGIME 



In the e-regime, it is necessary to evaluate the pion zero-modes, — (0,0), contribution non-perturbatively. It is 
convenient to integrate zero-mode out of the theory, leaving an effective field theory in terms of the non-zero modes. 
In the mixed regime, only the pion zero-modes are removed |30| . while the zero-modes of the kaon and eta are treated 
perturbatively. This can be done by rewriting the E operator as 



U = exp 







(5) 



The subscript z denotes zero-mode operators, while the operators with a hat are operators whose contribution can be 
treated perturbatively in the e-expansion. When integrating over the zero-modes it is convenient to write the operator 
U in terms of hyperspherical coordinates. 



^cos(V') -I- i cos{9) sin(V') sin(6') sin(?/')e** 0^ 



U 



\ 



sin(0) sin(V')e cos(V') — « cos(6') sin(?/') 

ly 



(6) 



When constructing the Lagrangian that is invariant under chiral transformations, it is advantageous to define the 
operator ^ = Under the redefinition of Eq. i5\, one finds, 



(7) 



where the definitions = iuWutU and V = VutU Wi^ have bee n implicitly introduced. When integrating 
over U, one may substitute Af^ = VAf_iV^ and = VVf_,V^ + iVd^V^ . The results presented here will be truncated 
at ©(e^), in which case one can safely make the following approximation: 



Mdf'M 



+ Oie'). (8) 



The only contribution to the heavy meson mass that originates from the zero- modes integration appears at ©(e"*), 
and comes from the second hne in Eq. Q: 



8Cm = ^Tr [HaH, [iMS, + h.c\^] + ^Tr (^X? + h.c.\,] 

~ ^Tr[H,H, + (cos(2^)m2) + ^TriHsHs] {2ml " 



(72 



(9) 
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A. 



Tr[HaHa] {2ml + (-1 + 2cos(2^)) ml) 



In order to evaluate the contribution of this term to the heavy meson mass, the ^-dependence in this expression must 
be integrated out using the non-pcrturbative weight arising from the last term in Eq. ([2| . It is convenient to perform 
this integral by analytically continuing to Euclidean time t — >■ —it : 



J VU^ exp J d^x (JjTr [MU^ + h.c] + SCm 



VU^ exp 
X{s) exp 



4 I d^x ((ai + 2a,) {p:^P: + pt . p^) _ (p*tp^* + pt . p^ 



(10) 



where s = \pm1j3lj^ , and X(s) can be expressed in terms of the modified Bessel function /i(2s) of the first kind, 

X{s)= /^dVcos2(i^)sin2(V;)e2^™^(2^) = ^^^, ( ^d^ cos''{i^)sin\^) e^'''°<^'^\os{2i^) ^ . (11) 



Since M has no spin structure, it results in the same shift for both the masses of the pseudoscalar and vector fields, 
yet it explicitly breaks the SU(3)v symmetry, 
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^'(^^ ■ '''2{2ml-ml), ^M(„ p. = 2^^^^ + ^^2(2m?, - rn^). (12) 
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This analysis introduces a volume dependence to the mass of the non-zero pion modes, as well as for the K^s and ry's. 
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(14) 



" 3 ^ 3 ''GXis) 3' 
After performing the integration over the zero-modes, the finite-volume contribution from the remaining degrees of 
freedom can be evaluated perturbatively. The finite volume Feynman diagrams can be evaluated in the standard way, 
where the integral is replaced by a sum over discretized four momenta and the zero mode is explicitly excluded in the 



pion loops |18| . An outline of the methods used in performing these sums is discussed in the appendix. Eq. (12 1 also 
includes a volume independent SU(3)y symmetry breaking term, 



S'q EE {6Mp^^ - (5Mp)(i^^) = 2ai 



{2ml - ml) 
Av 



(15) 



6 

This results in a shift on the bare mass of the strange pseudoscalar meson, Sq, and the strange vector meson 

A^^Aq + <S^. (16) 

At leading order, one may associate 5q with physical value of the splitting between the isospin doublet P and strange 
pseudoscalar Pg, which is on the order of 100 MeV for both the charm and bottom mesons respectively. For the sake 
of generality, both 6q and Aq will assume the same power counting as Aq ^ C(e^)- 

" •«.. 

/ \ 
' 

FIG. 1: e"^ contribution to the pseudoscalar heavy meson mass. The solid line corresponds to the heavy 
pseudoscalar, the double line denotes a vector meson, and the dashed line represents a Goldstone boson. 



RESULTS 

In taking the isospin limit, the pseudoscalar pair P — {P^, P^} will receive the same mass contribution. P will 
denote the isospin pair and Pg will denote the strange-heavy meson. In order to formally categorize the different 
terms contributing to the mass, it is important to consider the ratio Aqcd/wq ^ ©(e"). The most relevant cases are 
the following: 

(i) : AqcD/mq ^ ©(e^) (ii) : AqcdMq - 0{e), (17) 

corresponding to the static limit and LO heavy quark mass corrections, respectively. For simplicity, the expressions 
below will include finite LO heavy quark mass corrections. The static limit can easily be obtained by taking the 
toq oo limit (note Aq oc TOq^ — > 0). The individual diagrams contributing to Mp are written in the appendix. 
The notation SMp = Mp{L) — Mp{oo) is used to denote the finite-volume dependence of the mass. 

1. SU(3) HMxPT 

The SU(3) volume dependence of the P and Pg masses up to and including 0{e^) is found by adding the finite- 
volume contributions from the self-energy diagrams depicted in Figs. ([T]) and ([2]), where the Goldstone bosons can be 
pions, kaons and etas. 
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(M) m (c) 



FIG. 2: contribution to the pseudoscalar heavy meson mass, (a) Denotes the zero-modes contribution, (b) 
Goldstone bosons loops originating from four-point vertices, (c) Incorporates operators that contribute the heavy 
flavor symmetry breaking corrections to the PP*7r vertex. 
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where m'^ = 4rrij^/3, the discrete sums ci, C4, A/i, and /C; are defined in Eq. ( |A.20| , ( |A.24[A.25[ ), and Eq. ( [T5][T6| ) have 
been used. Similarly, the C(e"^) and ©(e"*) corrections to the vector meson masses are depicted in Fig. ^ and Fig. 



(a) 



(b) 



FIG. 3: contribution to the heavy vector meson mass. 



El), respectively. In the static limit pseudoscalar and the vector mesons are degenerate, therefore it is only necessary 
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to evaluate the volume dependence of the hyperfine splitting: 

SMp, - Syip = (^^^ + (mlM,{mK,L) + ^-mlMdm,,, L)^^ (20) 

5^Aq / ci 4 7V2(mx,i) , 2 7V2(to^,L)\ 



2/2 I 2^L2 3 ^1/2^5/2 9 ^1/2^5/2 



K 
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FIG. 4: e"* contribution to vector meson mass. 



2. SU(2) HMxPT 

In SU(2) chiral perturbation theory, the kaons and eta decouple from the theory. Only integrals including pions 
depicted in Figs. Q, ([s]), and Q contribute to the volume dependence of the masses, 

3 /„ (2) , „(2)\ C4 



(2) (2) /„(2)\2a(2) 



Q Ci 



(23) 



^mQ/2L3 /2 4^L2' 

where an additional superscript has been introduced in order to explicitly distinguish the SU(2) LECs from those 
contributing to the SU(3) theory. 

Note, these results have been derived assuming Aq ^ 0{e^), which should be expected to be the case for the charm 
mesons. For the bottom mesons one should expect Ab ~ C>{e^). This would move finite volumes effects related to 
this coupling to 0{e^), displacing them outside the scope of this calculation. 
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FIG. 5: In black is the fake data with corresponding error bars of approximately 0.5%. The red line is the fit of Eq. 



(24 1 to the data, and the filled area denotes the error associated with the fit of the LECs and Moo- The volume 



chosen for the data depicted in Fig. (6a) is 2 fm, and the pion mass chosen for the data depicted in Fig. (6c) is 0.2 
GeV. 



ANALYSIS AND DISCUSSION 

The results presented in the previous section allow determination of LECs that play an important role in heavy-light 
meson scattering. In order to evaluate the LECs, one must fit the expressions Moo + SM{L, ra-^) to LQCD results of 
the heavy meson masses for different volumes and pion masses that fall within the e-regime, where Moo is the infinite 
volume mass of the heavy meson and (5M(L,r7i^) is the finite volume contribution described by Eqs. 

The corresponding LQCD calculation has not been performed yet. Nevertheless, the uncertainty of the LECs as a 
function of the standard deviation of the heavy meson masses can be estimated. In order to do this I have analyzed 
fake data for the heavy meson mass. A data set was generated that follows the trend predicted by Moo + (5M(L, rriT^); 
this required inputing randomly-generated LECs. Additional volume and pion mass dependent terms were added 
to Moo + (5M(L,m„) in order to simulate the 0{t^) corrections. The exact form of these terms is irrelevant for the 
discussion at hand. Each heavy meson mass has been given a corresponding uncertainty. Lastly, the set is fit to 
Moo + (5M(L, m7r) in order to reproduce the randomly-generated LECs. Due to the larger number of undetermined 
LECs and the larger expansion parameters for SU(3) HM^PT (mx/A^^, m^/A^), the following discussion will focus 
on SU(2) HMxPT. 

depends on three parameters: m^, mq, and L. In order to determine a = {a\ + 2a\ )/2A^ 
and 7 = (27}"^'' -t- 72^^)/Ap(,, which contribute to the NLO calculation of Ptt scattering lengths [HIIZIj it is convenient 
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FIG. 6: Prediction for the level of precision for determining {g, 7, cr} by fitting Eq. (24) to six heavy mesons masses, 
with an uncertainty ranging from 0.5% to 0.01%, and each located at a different (L, m^r) that fall within the 
e-regime. The LECs used in generating the fake data were allowed to vary Moo = {4.5, 5.5} GeV, / — 130 MeV, 
g = {0.5, 1.5}, 7 = {0.5, 1.5} GeV"\ and a = {0.5, 1.5} GeV^\ 



to consider the static limit, toq ~ 00. In this limit the volume dependence of the mass at ©(e^) simplifies to, 

3 



2 X'(s) 3 C4 



(24) 



Since the level of precision with which the LECs can be determined depends on their magnitude, the parameters 
were varied Moo = {4.5,5.5} GeV, g = {0.5,1.5}, 7 = {0.5,1.5} GeV"\ and a = {0.5,1.5} GeV"\ and the pion 
decay constant was fixed a,t f — 130 MeV. For each value of the parameters, a set of six heavy meson masses was 
generated. Each heavy meson mass was calculated with different values of the volume and the pion mass that fall 
within the e-regime, and each heavy meson mass was given an uncertainty ranging from 0.5% to 0.01%. Fig. ([s]) 
serves as an example of one of these sets. The set depicted in Fig. ^ was generated using the center values for 
the parameters and each heavy meson mass was given an uncertainty of approximately 0.5%. The estimate of the 
expected fractional standard deviation of the LECs as a function of the uncertainty of the masses is plotted in Fig. 



11 

^ as the shaded region. The range of possible standard deviation for the LECs for a given uncertainty of the heavy 
meson mass manifests the fact that the precision with which these LECs can be determined depend on their absolute 
value. 

CONCLUSION 

HMxPT is the EFT for calculating strong-interaction quantities of heavy mesons. Currently, HMxPT is limited by 
the poor determination of the LECs of the theory. In particular, the LECs in the lagrangian discussed in this paper, 
Eqs. ([2]&|4]), contribute to the evaluation of scattering lengths, and are currently known to within a factor of three 
0[2!- 

These LECs can be evaluated from LQCD calculations. One way to achieve this is to extract the LECS from the 
volume dependance of the heavy mesons masses, since these finite volume effects are parametrized by the LECs of 
the theory. With this in mind, the finite-volume dependence of the heavy pseudoscalar and vector meson masses in 
the e-regime of HMxPT have been calculated to 0{e'^). In the e-regime, LQCD calculations can be performed at the 
physical point of QCD (m^ « 140 MeV) if volumes are small (L < 4fm). 

Lastly, it was shown that with a set of six heavy mesons masses extracted from LQCD with 0.1% uncertainties, 
the SU(2) NLO LECs a = (af ' + 2cr^^^)/2A^ and 7 = (27^^ + 7^^^)/A^ can be determined within the 10% level of 
precision. This level of precision shows the power of this method. In order to calculate the heavy-light scattering 
lengths, it is also necessary to determine the linear combination 274^'' + 73^'' [6l[7], which could be determined from 
the volume dependance of the heavy meson mass at 0{e^). 
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APPENDIX 



Finite- volume Feynman diagrams can be performed by replacing integrals with sums over discretized four momenta 
[18]. The sums contributing to the calculation of the heavy meson mass to 0(6^) are 



A(mi,A,L,/3) = — ^ y ^ ^Ar4 -, (A.l) 

/ 27rno -.2 

B(mi,L,/3) = — ^ V -Ar^ ^, (A.2) 



1 / 27ra\2 

C(mi,L,/3) = — ^ V ^At4 -. (A.3) 

where uj is the external energy, and mi denotes the light meson mass (m^, mj^, m^). In the mixed regime, the tt 
and {K, ij} loops must be treated separately. Due to the field convention, the corrections to the mass are defined as 
|n(a; = 0) + ^a<^n(a; = 0), where H is the sum of the amputated self-energy diagrams, and dp is the bare residual 
mass of the heavy meson. The superscripts of the terms below denote the order at which they contribute in the 
e-expansion. The ©(e'^) correction to the pseudoscalar mass, depicted in Fig. ([l]), is 



1 ^25^3^ ^ 1 

g^A(m., A, L, 13) = ^A(0, A, L, (3) + ©(e^), (A.4) 
M(3)(mi^,m,„A^L,/3) = ^A(mK, A^ L, /3) + (m^, A^ i, ^) (A.5) 



where 5 = 5-1- The first 0(6 ) contribution comes from integrating out the zero-modes using Eq. (12 1, depicted 



by Fig. (ga)), 

MW(m.,L,/3) = ^ia,+2a,)^. (A.6) 

The second graph. Fig. ([2|^b)), corresponds to the four-point vertex contribution to the mass, 

(m.,i,/3) = -5|±^B(0,L,/?) 
Mf^' {mK,m,,L, P) = -l^^B(mK, L, /3) - ^^^^^ (B(mK, L, /3) + C(m;^, L, /?)) 
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The third diagram, Fig. ([2][c)), comes from the three-point vertex corrections in the Lagrangian, and it results in the 
following contribution to the P meson mass: 

MW(m,, A, L, p) = -6-^A(0, A, L, p) + 0{e'), 

M^^\mK,m^A\L,(3) = --i|^A(mK, A^ L, /3) - x^f^A (^mK,A\L,(3] . (A.8) 

In order to evaluate the temporal sum, the Abel-Plana formula will be used: 

1 ■^-^ 27rn dz fiz) f(z) 

-al^f(-ff-'^ ^/(^:) - 'Res( ^,^^^ _ -^ jllowcrplanc + ^Res( ^_^^^ _ ^ ) [uppcrplane- (A.9) 

The spatial sum can be performed using Poisson's Resummation formula, 

n ^ ^ ^ n n 

'- V ' 

^ ^ ^ ^ > Y. - — • (A.io) 



As carefully discussed in Ref. |34| . the finite temperature contributions in the e- regime are usually heavily sup- 
pressed, and such is the case in all the integrals considered here. This allows one to safely neglect finite temperature 



terms. With this, one can extract the volume dependence of the sum in Eq. ( A.l I as 



5A(mi, A, L, /?) = A(mi, A, L, /3) — A(mi, A, i — > oo, /3 — > oo) 
1. e- Regime Integrals: mi = 



In the case that the sum arises from a pion loop, one can take the chiral limit and substitute \/k^ + mf — > kg in 

the above expression. Corrections to this approximation will result in 0{e^) contributions to the heavy meson masses. 

All the integrals can be preformed using the following generating formula: 

T(a, A) . r^fco^e-- = Ci(«A)^ ^ cos(aA) (| Si(.A)) ^ 
7q Kq + A A A 

poo h.2m a2m 

I"(a,A)^y^ dfc„_^e-^-o.^_I(,,A). (A.13) 



Where Ci(a;) = 7 + log(x) + "^ dt and Si(x) = rft are the geometric integrals. From Eq. (A.13), it 

follows: 

(5A(0,A,L,/3) =^ V ^— ( ~ Ci(nLA)sm(nLA) + cos(nLA) (sUnLA) - -)Y (A.14) 
47r^ nLA \ nLA V 2 / / 
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Assuming AL e, it is possible to expand about Ai = 0. In this limit the sum may be approximated as an integral 
over the variable z = nLA, [35], 

L^A^ H ~T\ \ ~T1\ ^ Ci(?^LA) sin(?iLA) + cos(nLA) (si(nLA) - J) ) 

Ci(^) sin(z) + cos(z) (si(z) - 



An 



27r" 



(A.15) 



At leading order this matches to the approximation made in Ref. ^35j for the same integral. The next term in the 
expansion comes from taking a derivative with respect to a' = LA 



d_ 



-a'^ ^— { — Ciina') sin(na') + cos(na') f Si(na') — - 

' ^-^ na' \ na' \ 2 



y- 



2n 



Expanding about LA leads to a 0{e ) approximation of Eq. (A. 14 1 



2L3 ^ 8n7rL2 



2L^ SttL^ 

0(e3) 0(€-*) 



where Eq. (A. 20 1 has been used. To 0{e ) the remaining integrals are 



47rL 



E 



27r 



C4 



2^2L* ^ n 



27r2L-^ 



5C(m.,L,/3) = --^5:ij+0(e 



27r2L^ 



(A.16) 



(A.17) 



+ 0(6*^), (A.18) 



(A.19) 



In writing out the full expression of the masses, it is important to note that 6M{0,L,/3) + (5C(0,L,/3) = + 0{e^). 
Two previously calculated sums have been used 



Cl 



y- 



-2.8372974 



C4 



y- 



16.532315. 



(A.20) 



2. p- Regime Integrals: m\ = {mKjm^} 



In the p-regime, the light meson mass is comparable to the lowest non-zero momentum mi/A^ ^ 2tt/LA^ ^ e. In 
this regime the small mass approximations used in the previous sections are no longer valid. One must perform the 



integral in Eq. (A. 11) without taking the chiral limit. Although this integral cannot be evaluated exactly, in the 



A — limit the integral is dominated by small values of ko. In this case, the argument in the exponential can be 
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approximated as ^/kQ +mf = mi + 



2mi Sm? 



SA{mi,A,L,l3) = A X! 

e ' 



dko 1 kg + nil ^ 



1 —nmiL 2-n 



2tt kl + A2 47rL 



nLkQ 



'"-1 \ ' 



A 



3 + 9minL — [m\nLY 
N2{mi,L)+0{e>), 



where the definition in Eq. (A.24) have been used. The remaining integrals can be performed exactly, 



(mi,i,/3) 



1 



E 



dki 



o,.2e 



2tt 



^0" 



167r2L 



^ (ii:3(TOi7iL) - KiimynL)) + ©(e^), 



TO, 



^/Ci(to:L)+0(.^) 

,3 

' ^ (KsiminL) + 3Ki{minL)) + ©(e^) 



TO, 



167r2L 



167r2L 



In writing Eqs. ( A.21|A!23 ) the following dimensionless functions were used, 



A/'i(toi,L) = 



-nnii L 



7V2(mi,L) = ^ 



3 + 9m,inL — {m,\nLY 
647r2n5/2 



K,i{m\L) — {K^{m\nL) — Ki{'m\nL)) IC-zimiL) — —4 Ki{m\nL), 



(A.21) 



(A.22) 



(/Ci(miL)-/C2(TOiL)). (A.23) 



2^g-nm,L (A.24) 

(A.25) 



where Ka are the modified Bessel functions of the second kind. Finally, by adding the contributions from Eqs. 
(A.4|A.8) and substituting the expressions of the respective sums, one arrives at Eq. (18). 
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